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We derive some new constraints on single-field inflation from the Wilkinson Microwave Anisotropy 
Probe 3-year data combined with the Sloan Luminous Red Galaxy survey. Our work differs from 
previous analyses by focusing only on the observable part of the inflaton potential, or in other words, 
by making absolutely no assumption about extrapolation of the potential from its observable region 
to its minimum (i.e., about the branch of the potential responsible for the last ~50 inflationary 
e- folds). We only assume that inflation starts at least a few e- folds before the observable Universe 
leaves the Hubble radius, and that the inflaton rolls down a monotonic and regular potential, with 
no sharp features or phase transitions. We Taylor-expand the inflaton potential at order v = 2, 3 or 
4 in the vicinity of the pivot scale, compute the primordial spectra of scalar and tensor perturbations 
numerically and fit the data. For v > 2, a large fraction of the allowed models is found to produce 
a large negative running of the scalar tilt, and to fall in a region of parameter space where the 
second-order slow-roll formalism is strongly inaccurate. We release a code for the computation of 
inflationary perturbations which is compatible with COSMOMC. 

PACS numbers: 98.80.Cq 



Cosmological inflation is known to be a successful 
paradigm providing self-consistent initial conditions to 
the standard cosmological scenario [E H, [E 0, [E 0] and 
explaining the generation of primordial cosmological per- 
turbations @, IE 0, EE EE E3 EE G3 ■ The distribution 
of Cosmic Microwave Background (CMB) anisotropics, 
as observed for instance by the Wilkinson Microwave 
Anisotropy Probe (WMAP) [IE EE EE El, is compati- 
ble with the simplest class of inflationary models called 
single-field inflation. 

The definition of single-field inflation is not unique: 
for instance, some authors consider hybrid inflation 
[IE HE EH as a multi-field model, since it involves one 
scalar field in addition to the inflaton field (the role of 
the second field being to trigger the end of inflation). 
In this work, we call single-field inflation any model in 
which the observable primordial spectrum of scalar and 
tensor metric perturbations can be computed using the 
equation of motion of a single field. This definition does 
include usual models of hybrid inflation. 

The goal of this paper is to derive from up-to-date cos- 
mological data some constraints on the scalar potential 
V(4>) of single-field inflation. This question has already 
been addressed in many interesting works since the publi- 
cation of WMAP 3-year results [jj, M HE EI HE HE EE 
HE HI, HE EE HI (see also [H for earlier results). Our 
approach is however different, since all these references 
assume either that the slow-roll formalism can be applied 
(at first or second order), or that the scalar potential can 
be extrapolated from the region directly constrained by 
the data till the end of inflation. We want to relax these 
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two restrictions simultaneously, and to derive constraints 
on the observable part of the inflaton potential under the 
only assumption that V(</>) is smooth enough for being 
Taylor-expanded at some low order in the region of in- 
terest. In this respect, our work is still not completely 
general and does not explore possible sharp features in 
the inflaton potential (see e.g. [IE EE] for recent propos- 
als). Throughout the abundant literature on the inflaton 
potential reconstruction, the work following the closest 
methodology to ours is the pre- WMAP paper of Grivell 
and Liddle [34| . 

The question of whether the slow-roll formalism can 
be safely employed or not is intimately related to the 
magnitude of a possible running of the tilt in the pri- 
mordial spectrum of curvature perturbations. In order 
to clarify this point, lets us first recall that the slow-roll 
formalism [2E EE E3 consists in employing analytical 
expressions for the primordial spectrum of curvature per- 
turbations Vn(k) and gravitational waves Vh(k). Such 
expressions hold in the limit in which the first and second 
logarithmic derivative of the Hubble parameter H with 
respect to the e-fold number N = In a remain smaller 
than one throughout the AN ~ 10 observable e-folds of 
inflation (i.e., over the period during which observable 
Fourier modes cross the Hubble radius). Deep inside this 
limit, the primordial spectra are given by 

V ^ k ) ~ 2 uTuum > ?V fc ) * , (1) 

irmp(dmH/dN) nm z p 

where the right-hand sides are evaluated at Hubble cross- 
ing. The first-order expression of the scalar/tensor tilts 
ns.T and tilt runnings as,T can be easily obtained by 
taking the derivative of the above expressions, using the 
slow-roll approximation d/dlnfc ~ d/dN. The derivation 
of higher-order ex pres sions is more involved (see e.g . |3 8l 

EE EE HE M MM M M M M M EE HE Ha M) ■ 
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FIG. 1: Observable region of the inflaton potential allowed at the 95% C.L. by the WMAP 3-year and Sloan Luminous Red 
Galaxy survey (SDSS-LRG) data, for a Taylor expansion of the potential at order v = 2 (top), v = 3 (middle) or v = 4 
(bottom), in the vicinity of the pivot value <j>*. In all diagrams the potentials are normalized to their value at the pivot scale 
4>*. For clarity, in the right diagrams, the field is expressed in units of l^/V^'l instead of mp, so all curves have by definition 
the same slope in cf> = In practice, these plots show the superposition of 95% of the potentials from our MCMC chains with 
the best likelihood (after removal of the burn-in phase). Each potential is plotted in a range [4>i,4>2] corresponding to Hubble 
exit for modes in the range [fci,fc2] = [2 x 1CP 4 ,0.1] Mpc -1 which is most constrained by the data. This corresponds roughly to 
a history of 6.2 e-folds. We only show here potentials with a negative slope, but their image under the cj> — > —cj> symmetry are 
equivalent solutions. At first sight, on the top left diagram, v = 2 potentials seem to have a non-zero third derivative, but this 
impression comes only from the superposition of many lines with varying length and slope. 



Current data clearly indicate that around the pivot 
scale at which the amplitudes, tilts and runnings are de- 
fined (usually, the median scale probed by the data), the 
tensor-to-scalar ratio is small and the scalar tilt is close 
to one. This is sufficient for proving that the two slow- 
roll conditions are well satisfied around the middle of 
the observable e-fold range. However, depending on the 
inflaton potential, higher derivatives d n \nH/dN n (with 



n > 3) could be large near the pivot scale, leading to a 
sizable tilt running as and eventually to a situation in 
which slow-roll would hold only marginally at the begin- 
ning and/or at the end of the observable e-fold range. 
This explains why the two issues of large \ots\ and slow- 
roll validity are closely related (as recently emphasized 
in H). 

Since models with a large running imply that the two 
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slow-roll conditions become nearly saturated near the 
ends of the observable potential range, a naive extrapola- 
tion would suggest that they cannot sustain inflation for 
much more than the observable 6N ~ 10 e-folds. How- 
ever, it is always conceptually easy to extrapolate the 
potential in order to get the necessary 50 or 60 infla- 
tionary e-folds after our observable universe has left the 
Hubble radius, or to obtain arbitrary long inflation be- 
fore that time. Potentials designed in that way might 
not have simple analytical expressions. This should not 
be a major concern e.g. for physicists trying to derive 
inflation from string theory, in which the landscape de- 
signed by the multi-dimensional scalar potential can be 
very complicated, leading a priori to any possible shape 
for the effective potential of the degree of freedom driving 
inflation. However, it is clear that models inducing large 
running are not as simple and minimalistic as those with 
negligible running. But since they are not excluded, they 
should still be considered in conservative works such as 
the present one. 

In section HI we will follow a conventional approach 
and fit directly the Taylor-expanded primordial spectra 
to the data. Like most other authors, we will conclude 
that: (i) the data provides absolutely no indication for 
as ^ 0, and (ii) given the current precision of the data, a 
large running is nevertheless still allowed. We will show 
that similar conclusions also apply to the running of the 
running (3 s- 

In section [H] which contains our main original results, 
we will fit directly the Taylor-expanded scalar potential 
of the inflaton to the data. Our reconstructed potentials 
are displayed in Fig. [TJ Unless we impose a "no-running 
theoretical prior" (i.e., the prejudice that inflation is deep 
inside the slow-roll regime) , our potentials will freely ex- 
plore the region in parameter space where the running 
(and eventually the running of the running) are as large 
as found in sectionQ] So, for self-consistency, we must for- 
get about the slow-roll formalism and compute the exact 
primordial spectra numerically (as Ref. [26j did for vari- 
ous specific exp ressions of the potentials). In a very nice 
work, Ref. [54[ gave a few examples of scalar potentials 
leading to the largest \as\ values allowed by the data, and 
showed that even in these cases the second-order slow-roll 
formalism, although inaccurate, remains a reasonable ap- 
proximation. In the present systematic analysis, which 
explores the full parameter space of smooth inflationary 
potentials, we will sec that this conclusion docs not apply 
in all cases. 



I. FITTING THE PRIMORDIAL SPECTRUM 

Primordial spectrum parametrization. The usual way 
of testing inflationary models without making too many 
assumptions on the inflaton potential is to fit some 
smooth scalar/tensor primordial spectra, parametrized 



as a Taylor expansion of ln'P with respect to lnfc, 



In 



Vn(k) 



fc as 2 & 
[n s - 1) In — + — In — 



0S, 3 



k 



(2) 

and the same holds for Vh as a function of tit, cut and 
Pt- In single-field inflation, the coefficients of the scalar 
and tensor spectra are related through the approximate 
self-consistency relation 



dlnV h (k) _ 1 V h (k) 
d\nk ~ 8Vn(k) 



(3) 



which follows trivially from Eq. ([T]) and becomes exact 
deep in the slow- roll limit. The sensitivity of current data 
to gravitational waves is very low, with loose constraints 
on the shape of Vh- So, even if the slow-roll formalism 
might become inaccurate in some cases, the data can be 
fitted assuming that Eq. j3]) is exact. In other words, 
for practical purposes, we can safely use the hierarchy of 
relations derived from Eq.([3]), 



TIT 



-r/8, 



a T = n T [n.T - n s + 1], 



etc. 



(4) 



where r = Vh{k*)/Vn(k*)- So, if we decide to Taylor 
expand the scalar spectrum with p independent coeffi- 
cients, the total number of free inflationary parameters 
in the problem is p+1, including the tensor-to-scalar am- 
plitude ratio at the pivot scale. 

In principle, p should be chosen according to Occam's 
razor: when increasing p does not improve sufficiently the 
goodncss-of-fit, it is time to stop. In a Bayesian analy- 
sis, this question is addressed by the computation of the 
Bayesian evidence [H, H3, [H, [H, Ififj, EE IH3|- However, 
when the evidence does not vary significantly as a func- 
tion of p, the decision of stopping the expansion remains 
a personal choice to some extent, and more conservative 
works should consider higher p values. 

The issue of varying p is important in two respects: 
first, one needs to know how many independent informa- 
tions the data is providing, i.e., how smooth/complicated 
the inflaton potential needs to be (for addressing this is- 
sue, one could also perform a principal component analy- 
sis [63|); second, it is useful to know whether the bounds 
on a given cosmological/inflationary parameter 0j are in- 
dependent of p, or subject to variations when p increases, 
due to the appearance of new parameter degeneracies. 

Results. In order to address these two points, we per- 
formed some global parameter fits using the public code 
COSMOMC [6J|, with p varying from two (scalar ampli- 
tude and tilt) to four (including the tilt running, as well 
as the running of the running). Our results are sum- 
marized in Tabic HI The relative Bayesian evidence of 
each model can be easily computed, since the models 
are nested inside each other [53]. However, this calcu- 
lation forces us to choose some explicit Bayesian priors 
for as and (3s- Pushing inflation to its limits, we no- 
tice that the second slow-roll parameter can in princi- 
ple vary between plus and minus one, so the scalar tilt 
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FIG. 2: Probability distribution of cosmological and inflationary parameters for the models of section U p — 2 (green/light 
solid), p — 3 (black dashed), and the models of section HT1 v = 2 (magenta dot-dashed), v = 3 (blue dotted), v = 4 (red/dark 
solid). For the runs of section U the free parameters (with flat priors) are the first eight parameters; the corresponding 
probability for the potential parameters are derived from second-order slow-roll formulae (involving V* to Vl" , so the inferred 
value of V"" remains undetermined). Instead, for the runs of section [TTJ the free parameters are the first five and the last four; 



the amplitude parameter in the fifth plot is then defined as In 



10 
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the corresponding r, ns and as are derived from 
the exact primordial spectra. The data consists of the WMAP 3-year results [lH [l6l. IvH, [lSj and the SDSS LRG spectrum [55( [. 



:SV,'- 



could take any value between zero and two. Extreme 
runnings could be observed in ad-hoc inflationary mod- 
els such that during the observable e-folds, correspond- 
ing to four decades in k space, ns evolves from to 2 
or vice- versa. So, the as prior can be chosen to be a 
top-hat centered on zero with Acts = 4/ln(10 4 ) ~ 0.4. 
Similarly, extreme values of Ps correspond to ns passing 
through the sequence 0-2-0 or 2-0-2. This leads to a prior 
width A/? s = 16/ln 2 (10 4 ) ~ 0.2. With such priors, the 
Bayesian evidence E increases by a factor 

§^ = [V{a s = 0) Ac,]- 1 = [2.1 x 0.4]- 1 = 1.2 (5) 

when as is added, and again by 

Ep± = [<p(j3 s = o) A/3,]- 1 = [6.2 x 0.2]" 1 - 1.2 (6) 



when (3s is introduced. These numbers are too close to 
one for drawing definite conclusions: the extra parame- 
ters are neither required, neither disfavored by Occam's 
Razor. 

TableQ]shows that adding as has a small impact on the 
probability distribution of flbh 2 , r, r, Vn(k*) and ns, as 
found in previous works. However, it is reassuring to note 
that adding j3s leave all bounds perfectly stable, except 
for a small shift to higher n s values. This suggests that 
including a few higher derivatives beyond as does not 
open new parameter degeneracies (this conclusion would 
probably break if the number of free parameters becomes 
much larger). Figs. [2 [3] show the likelihood distribution 
of each parameter as well as some two-dimensional con- 
fidence regions for models p = 2 (green lines) and p = 3 
(black lines). 
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FIG. 3: Two-dimensional 68% and 95% confidence level con- 
tours based on WMAP 3-year and the SDSS LRG spectrum 
for the parameters describing the primordial spectra, obtained 
directly from the MCMC in the case of models p — 2 (green) 
and p = 3 (black) , or derived form the exact spectra for mod- 
els ii = 2 (magenta), v = 3 (blue), v = 4 (red). 



Parameter 


p = 2 


P = 3 


p — 4 




0.0226 ± 0.001 


0.022 ±0.001 


0.021 ±0.001 


O h 2 


0.109 ±0.004 


0.109 ±0.004 


0.109 ±0.005 


e 


1.041 ± 0.004 


1.040 ± 0.004 


1.041 ±0.004 


T 


0.08 ±0.01 


0.10 ±0.02 


0.10 ±0.02 


ln[10 10 P^] 


3.08 ± 0.06 


3.13 ±0.07 


3.13 ±0.07 


r 


< 0.13 


< 0.3 


< 0.3 


ns 


0.97 ±0.02 


1.00 ±0.04 


1.03 ±0.05 


as 





-0.07 ±0.04 


-0.07 ±0.04 










-0.04 ±0.04 


111 ^max 


2688.3 


2687.1 


2686.5 


E 


1 


1.2 


1.4 



TABLE I: Bayesian 68% confidence limits for ACDM infla- 
tionary models with p = 2, 3, 4 coefficients in the logarithmic 
Taylor expansion of the scalar primordial spectrum. The last 
lines show the maximum likelihood value and the Bayesian 
Evidence (relative to that of p = 2). The data consists in 
the WMAP 3-year results [H, El E3, El and the SDSS-LRG 
spectrum [55j], as implemented in COSMOMC [64| . 



Impact of small CMB multipoles. The smallest mul- 
tipolcs in the CMB temperature and polarization maps 
arc still controversial. In WMAP data (as well as in 
previous COBE data) , the temperature quadrupoles and 
octopoles are surprisingly small, while their orientations 
seem to be correlated (between each other and with the 
ecliptic plane). Many authors have been investigating 
possible foregrounds or systcmatics which could affect 



these small multipoles (see e.g. [63, ISs |67| and refer- 
ences therein). So, it is legitimate to study whether the 
quadrupole and octopole data have a significant impact 
on our bounds for the primordial spectrum parameters 
(a priori, these low temperature multipoles could be par- 
tially responsible for the preferred negative value of the 
tilt). 

We repeated the p = 3 analysis after cutting the tem- 
perature and polarization data at I = 2, 3. We found that 
all probabilities arc essentially unchanged, including that 
for running (the mean value only moves from —0.67 to 
—0.66, which is not significant given the precision of the 
runs). We conclude that our results are independent of 
the robustness of low mutipole data. fn\ 

Impact of extra CMB data. There are also discussions 
about a possible small mismatch in the amplitude of the 
third acoustic peak probed on the one hand by WMAP, 
and on the other hand by Boomerang [68j or other small- 
scale CMB experiments. We repeated the p = 2 and 
p = 3 analysis including extra data from Boomerang [68| , 
ACBAR [69| and CBI 0. The impact on inflationary 
parameter is found to be very small, although in the p = 3 
case the bound on r gets weaker by 20% and the preferred 
value of as goes down by the same amount (i.e., the case 
for negative running becomes slightly stronger). Other 
bounds are essentially unchanged. In what follows, we 
will not include these data sets anymore. 

Expectations for the infiaton potential. In the next 
section, we will directly fit the infiaton potential V (</>), 
parametrized as a Taylor expansion near the value 
corresponding to Hubble crossing for the pivot scale fc*. 
We expect that a global fit with V(cf>) expanded at order v 
will provide the same qualitative features as the previous 
power spectrum fit of order p = v: 

• order v — 2 (V'" = 0) should be sufficient for ex- 
plaining the data, and will not lead to significant 
running or deviation from slow-roll. Indeed, the 
smallness of r and |ng — 1| guarantees that the two 
slow-roll conditions are well satisfied at least near 

In addition, with V'" = 0, they should remain 
well satisfied on the edges of the spectrum, and no 
significant running can be generated. 

• order v > 3 (V"" = 0) should not be required by 
the data, but remains interesting since it will ex- 
plore the possibility of large running and shift the 
other parameter distributions as in the previous 
p = 3 case. The slow-roll parameters could then 
become large near the edges, so it is necessary to 
compute the spectra numerically rather than using 
any slow-roll approximation. 

The results of the next section will confirm this expecta- 
tion, and prove that order v = 4 is necessary for exploring 
the full range of a s probed by the p = 3 run. 
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II. FITTING THE SCALAR POTENTIAL 

Computing the power spectra numerically. In order to 
fit directly the inflaton potential, we wrote a new COS- 
momc module which computes the scalar and tensor pri- 
mordial spectra exactly, for any given function V(0— 0*). 
This module can be downloaded from the website 
http : //wwwlapp . in2p3 . f r/~lesgourgues/inf lation/ 
and easily implemented into COSMOMC. 

In its present form, our code is not designed for mod- 
els with very strong deviations from slow-roll. For such 
extreme models, a given function V(4> — 0*) would not 
lead to a unique set of primordial spectra Viz(k), Vh(k): 
the result would depend on the initial conditions in phase 
space. We decide to limit ourselves to models such that 
throughout the observable range, the field remains close 
to the attractor solution for which is a unique function 
0. In this case, a given function V(<fr — 0*) does lead 
to unique primordial spectra, and we do not need to in- 
troduce an extra parameter 0; n i. Since the goal of this 
paper is to test inflationary potentials leading to smooth 
primordial spectra, this restriction is sufficient. In par- 
ticular, it enables to explore models for which the run- 
ning as is large, deviations from slow-roll are significant, 
and analytical derivations of the spectra are inaccurate. 
However, our code cannot deal with the case in which 
inflation starts just when our observable universe crosses 
the horizon (for which 0j n i would be a crucial extra free 
parameter). 

In COSMOMC, we fix once and for all the value of the 
pivot scale fc» = 0.01 Mpc -1 . Then, for each function 
V((j) — 0*) passed by COSMOMC, our code computes the 
spectra Vn{k), Vh{k) within the range [fc m i n , fc max ] = 
[5 x 10~ 6 ,5] Mpc -1 needed by CAMB, imposing that 
aH = fc* when = 0*. So, the code finds the attractor 
solution around = 0*, computes -ff* and normalizes 
the scale factor so that a* = fc*/-ff*. Then, each mode 
is integrated numerically for k/aH varying between two 
adjustable ratios: here, 50 and 1/50. So, the earliest 
(resp. latest) time considered in the code is that when 
kmm/o-H = 50 (resp. k max /aH = 1/50), which in the 
attractor solution uniquely determines extreme values of 
(0 — 0*) according to some potential. In the code this is 
translated to demanding that aH grows according to the 
aforementioned ratios: by 50fc*/fc m i n before = 0*, and 
by 50fc max /fc* afterwards. Hence, one of the preliminary 
tasks of the code is to find the earliest time. If by then, 
a unique attractor solution for the background field can- 
not be found within a given accuracy (10% for 0i n i), the 
model is rejected. So, we implicitly assume that inflation 
starts at least a few e-folds before the present Hubble 
scale exits the horizon. In addition, we impose a posi- 
tive, monotonic potential and an accelerating scale factor 
during the period of interest. This prescription discards 
any models with a bump in the inflaton potential or a 
short disruption of inflation, that could produce sharp 
features in the power spectra. 

As a result of the chosen method, the potential is 



slightly extrapolated beyond the observable window, in 
order to reach the mentioned conditions for the begin- 
ning and ending of the numerical integration. Although 
this seems to be in contradiction with the purpose of this 
paper, i.e. to probe only the observable potential, this 
extrapolation cannot be avoided if we want to keep the 
number of free parameters as small as possible. Note that 
the range of extrapolation is still very small in compari- 
son with an extrapolation over the full duration of infla- 
tion after the observable modes have exited the Hubble 
radius. 

In this approach we need not make any assumption 
about reheating and the duration of the radiation era. As 
explained in Ref . [3l| , the evolution during reheating de- 
termines the redshift z* at which presently observed per- 
turbations left the Hubble radius during inflation. Prob- 
ing only the observable window of perturbations, we are 
allowed to let the subsequent evolution of inflation and 
reheating, hence the number of c-folds and thereby the 
redshift z*, be unknown. 

Parametrization. The inflaton potential is Taylor ex- 
panded up to a fixed order, and we let COSMOMC probe 
different values for the derivatives of the inflaton po- 
tential at the pivot scale. Since Monte Carlo Markov 
Chains (MCMC) converge faster if the probed parame- 
ters are nearly Gaussian distributed, in fact we recom- 
bine to potential parameters in such a way as to probe 
nearly Gaussian combinations. These combinations are 
inspired by the slow-roll expression of the spectral pa- 
rameters (T > -R,(k*), ns, ots and r) as a function of the 

potential. We use an amplitude parameter v X* m 6 , 

which is equal to Viz^k^) at leading order in a slow-roll 
expansion. The other spectral parameters consist of lin- 
ear combinations of (K'/V*) 2 , V?/V, (K"7K)(K7K) 
and (Vf /V*){Vi/V*f. Hence it is most likely to find 
nearly Gaussian shapes for these products in stead of the 
sole potential derivatives. For the actual expressions for 
the spectral parameters in terms of the inflaton potential, 
we refer the reader e.g. to section IV of (44[. 

In order to compare the results for the runs with v = 
2,3,4 with those of the previous section, we calculate 
the spectral parameters of each model numerically. The 
other way around, we also invert the slow-roll expansion 
in order to compare the v = 2, 3, 4 and p = 2, 3 models in 
potential-derivative space. Defining eo = H(Nj)/ 'H(N) 

and e n +i — d '^jy" ^ ' wncrc N = In ^- and Hi is some 
initial value of the Hubble factor, the inversion is given 
by 

r C\ (r 2 , ^ \ 

ei = 1 \-(n s -l)r\ 

16 16 V 8 J 

+o(r 3 ,(n s -l) 3 ,«|) , (7) 
e 2 = -(ns-^ + das-^-^ins-^ld-^) 
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Parameter 



In 



1287rlO iu V? 



V" 2 

-vr m p 

vt" Vi 4 
-^-■^r-mp 

v'J" (vl\ 

v * \ V 'J 



In £„ 



tion 



0.022 ± 0.001 0.022 ± 0.001 0.022 ± 0.001 

0.109 ± 0.004 0.109 ± 0.004 0.109 ± 0.004 

1.041 ± 0.004 1.041 ± 0.004 1.040 ± 0.004 

0.08 ±0.03 0.09 ±0.03 0.10 ±0.03 

3.06 ±0.06 3.07 ±0.06 3.11 ±0.08 



< 0.4 

0.1 ±0.5 






< 0.4 

-0.2 ±0.' 
8±5 




< 0.8 

0.4 ±0.9 
13 ±11 
200 ± 150 



2688.3 



2687.2 



2687.2 



TABLE II: Bayesian 68% confidence limits for ACDM infla- 
tionary models with a Taylor expansion of the inflaton po- 
tential at order v — 2, 3, 4 (with the primordial spectra com- 
puted numerically). The last line shows the maximum like- 
lihood value. The data consists of the WMAP 3-year re- 
sults EE BS E3, and the SDSS LRG s P ectrum HI' as 

implemented in COSMOMC [64| . 



-j (C 1 -l)±0(r a ,(n s -l) 3 ,4), 
(r 2 



± (ns - l)r — Sa s 
+0 (r\{n s ~lf ,a s ) , 



(9) 



where C x = ^ E + In 2 - 2 ~ -0.7296. The value of the 
potential and its derivatives can be expressed exactly in 
terms of the slow-roll parameters, which are listed up to 
the second derivative in 4J|. The third derivative reads 



v ,„ = V2mlH\ 



2e\- 



3e 2 £i 



The fourth derivative of the inflaton potential would be 
of a higher order in the slow-roll expansion. 

Results. The allowed ranges, parameter likelihoods 
and two-dimensional contours from all our runs are sum- 
marized respectively in Table HH Fig. [5] and Figs. [31 2] 
The allowed shape of primordial scalar and tensor spec- 
tra is shown in Figs. [5) 

First we ran a chain for the model at order v = 2. 
As expected, the results confirm those obtained fitting 
the spectral parameters up to order p = 2, which can be 
seen in Figs. [5] and 2] and the upper left chart in Fig. [31 
by comparing the magenta and green lines. This can be 
translated into the statement that fixing the running of 
the tilt to zero is almost equivalent to fixing third and 
higher derivatives of the potential to zero. The resulting 
bounds can be read from Table [Til an d the correlation 
between V£/V* and V" /V* is well accounted by the rela- 
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= 0.6 ±0.2 



(68%C.L. 



(11) 

Note that the numerically calculated running in the mod- 
els with v = 2 is not strictly zero but allows a very small 
region of nonzero running. Similarly, the derived bounds 
from the models with p = 2 on the potential parameters 
allow for very small regions of nonzero second and third 
derivative. This merely reflects the expansions in differ- 
ent paramctcrizations than an indication for running (or 
nonzero higher derivatives of the potential). 

Including a third derivative does allow for models to 
have a more significant running, which is clearly visible 
Fig. [5j Yet, as seen in Figs. [21 [3] and HI the models with 
v = 3 (blue line) do not explore the full range of param- 
eters which is indicated by the models with p = 3 (black 
line), in particular for as, and do not show as much a 
sign of degeneracy between V"'Vl/V 2 and VI 2 as the de- 
rived potential derivatives from the models with p = 3, in 
Fig. [H The relation between V" and V' 2 remains almost 
unchanged. Inversely, in Fig.[3Jwe see the same effect in 
spectral parameter space. 

The remaining discrepancy between models v = 3 and 
p = 3 led us to including the fourth derivative of the 
inflaton potential as a free parameter, i.e. v = 4. The 
resulting power spectra, shown in Fig. show a larger 
negative running than in the v = 3 case, even with sig- 
nificant running of the running on the largest scales. In 
Fig. [2] we see that the model v = 4 (red line) does probe 
the same range of runnings of the tilt as allowed in the 
model p = 3. Looking at two-dimensional projections, 
we see that the p — 3 and v — 4 contours are closer to 
each other in spectral parameter space (Fig. [3]) than in 
potential parameter space (Fig. [4]): this reflects the inac- 
curacy of second-order slow-roll expressions, as explained 
in the next paragraph. In the model v = 4 the range for 
the lower derivatives of the potential is slightly larger 
than in the models with v < 4, which has its repercus- 
sions embodied in slight degeneracies between the fourth 
derivative and the lower derivatives. 

Note that all figures containing information on the 
fourth derivative of the potential contain only the model 
v — 4 (red line) and not those with p = 2 or p = 3, since 
in the slow-roll approximation one would need to go to 
third order in order to infer V"" from the primordial 
spectrum. 

Finally, it is worth pointing out that the results at 
all orders in both paramctcrizations still allow for a flat 
(Harrison-Zel'dovich) spectrum at the 95% C.L., or for a 
linear potential at the 68% C.L. 

Precision of the slow-roll approximation. In Fig. \E\ we 
show the discrepancy between the numerical results for 
the spectral parameters (top left: r, top right: ns, bot- 
tom left: as) and those obtained using the slow-roll ap- 
proximation up to third order in the derivatives of the in- 
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FIG. 4: Two-dimensional 68% and 95% confidence level contours based on WMAP 3-year and the SDSS LRG spectrum, for 
the parameters describing the inflaton potential, obtained directly from the MCMC in the case of models v — 2 (magenta), 
v — 3 (blue), v = 4 (red), or derived from second-order formulas for models p — 2 (green), p = 3 (black). 



fiaton potential (second order in slow-roll parameters for 
ns and as, third order for r). The numerically calculated 
parameters can in this context be treated as exact, since 
they do not involve any approximation (within first-order 
cosmological perturbation theory), and remain perfectly 
stable when we increase the precision parameters of our 
code. As r naturally comes out at one order higher in 
the slow- roll expansion than n$ and as, the top left plot 
(r) shows less discrepancy than do the plots for ns and 
cvs- However, for large r there is a clear deviation from 
slow-roll in the results with v — 4 (red), up to ~7%. In 
the plot for ns it is clearly seen that second order slow 
roll is only accurate up to ~ 3% for the run with v = 3 
(blue region). This discrepancy is important, since the 
data constrains ns with a standard deviation a ~ 2%. 
When a fourth derivative of the inflaton potential is in- 
cluded, second-order slow roll becomes really inaccurate, 
with a typical error of 10% on ns, while the running can 
be wrongfully estimated by as much as Aas =0.1, i.e. 
three standard deviations given the current data. In a 
future work, it would therefore be useful to compute the 
next-order contributions to the running analytically (the 



bottom right diagram shows the quasi-linear dependence 
of as on the combination V^'V't/V?). 



III. CONCLUSIONS 

In this work, we derived some constraints on the in- 
flaton potential from up-to-date CMB and LSS data. 
Our CMB data consists in the WMAP 3-year measure- 
ment of the temperature and polarization power spec- 
trum. We did include the first (controversial) multipoles, 
after checking in section U that they do not have a sig- 
nificant impact on the determination of the primordial 
spectrum tilt and running. Our analysis differs from pre- 
vious works for several reasons. First, we directly fit 
the parameters describing the inflaton potential, instead 
of constraining first the primordial spectra, and recon- 
structing the inflaton potential afterwards. Second, we 
Taylor-expand the inflaton potential in the vicinity of the 
pivot scale at a rather high order (up to v = 4), and sec 
that such a high order is important e.g. for exploring 
all the parameter space allowed by the data in terms of 
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FIG. 5: The primordial spectrum for scalar perturbations (left) and tensor perturbations (right) allowed at the 95% C.L. by 
WMAP 3-year and the SDSS LRG data, for a Taylor expansion of the inflaton potential at order v = 2 (magenta/ light), v = 3 
(blue/dark) or v — 4 (red/medium). In practice, this plot shows the superposition of 95% of the spectra from our MCMC 
chains with the best likelihood (after removal of the burn-in phase). All these spectra are computed numerically, rescaled to 
one at the pivot value k», and displayed in the range which is most constrained by our data set. Note that the shapes of the 
two spectra are related to each other: so, the tensor spectrum is constrained through that of the scalar spectrum, and not 
directly by the data, which does not have the required sensitivity. 
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FIG. 6: Precision test of the second-order slow-roll expansion for models allowed at the 95% C.L. by WMAP 3-year and the 
SDSS LRG data, based on a Taylor expansion of the inflaton potential at order v = 2 (magenta/ light), v = 3 (blue/dark) or 
v = 4 (red/medium). For each model, we plot the spectral parameter r (top left), ns (top right) and as (bottom left) computed 
at the pivot scale with two methods: by deriving the primordial spectrum computed numerically (horizontal axis), or with the 
second-order slow-roll formalism. We checked that the scattering of the point away from the y=x axis reflects inaccuracies in 
the second-order slow-roll formalism rather than in our code. The bottom right plot shows the inaccuracy in the running as a 
function of V""V'1/V, S for the v = 4 model. 



running of the scalar spectrum tilt. Third, wc compute 
the scalar and tensor spectra for each model numerically, 
and find that for the models considered here this is im- 
portant, since the spectra derived from the second-order 
slow-roll formalism are inaccurate by the same order as 



the observational constraints themselves. 

However, the most important peculiarity of this work 
is our choice to focus only on the observable region of 
the inflaton potential, not making any assumption on the 
shape of the potential between the observable region and 
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the minimum close to which inflaton stops after approx- 
imately 50 e- folds (depending on the scale of inflation). 
This choice has a crucial impact on the results. If we did 
extrapolate the inflaton potential over 50 e-folds, keep- 
ing the same order in the Taylor expansion, our models 
would be more severely constrained, since the require- 
ment of 50 extra e-folds would kill many of the allowed 
potentials presented here[72|. We are perfectly aware of 
this, and wish to point out that this is one of two points 
of view, which are both equally sensible. 

From one point of view, if one works under the preju- 
dice that the inflaton potential should not be too compli- 
cated, then it is extremely relevant to consider the global 
shape of the potential and to throw away all models which 
cannot sustain 60 inflationary e-folds. Many papers use 
this approach, using sometimes Monte Carlo methods in 
which the potential (or the Hubble flow H(N)) is Taylor 
expanded over the 60 e-folds at high order. 

From another point of view, if one wants to address 
the question of what is strictly allowed by the data, then 
even a high-order Taylor expansion of the full potential 
sounds unsatisfactory for modeling all its possible varia- 
tions during such a long history as 60 c-folds (especially if 
one keeps in mind that some other fields could then play a 
role: triggering a phase transition, inducing complicated 
shapes as in the string-inspired landscape scenarios, etc.). 
On the contrary, in this philosophy, one should only try 
to parametrize the inflaton potential in the range probed 
by cosmological data, i.e. around six or seven e-folds. 
This is what we did here, with a Taylor expansion up to 
fourth order. 

The two approaches lead, of course, to radically differ- 
ent conclusions. For instance, in the first method, one 
would conclude that during the observational e-folds the 
inflaton must be deep in the slow-roll regime, since it 
is necessary to sustain a number of e-folds which is an 
order of magnitude higher. The running would then be 
very constrained (24] ]. In the second method, it is not a 
problem to satisfy slow-roll only marginally on the edges 



of the observable range. Even if ei grows dangerously 
close to one when cluster scales exit the Hubble radius, 
the potential could become much flatter afterwards, and 
sustain any desired amount of inflation. 

Our main results for the inflaton potential reconstruc- 
tion are summarized in Figs. [TJ [21 SI an d Table [TTJ We 
also showed up to what extent the slow-roll formalism 
reveals to be inaccurate in the current context in Fig. [5] 
This motivates possible future works concerning the next- 
order slow-roll expressions. 

Following the same approach, this work could be im- 
proved by adding more large-scale structure data e.g. 
from Lyman-a forests or weak lensing, which have a good 
power for further constraining the primordial spectrum 
on smaller scales than the SDSS LRG data. Here we 
choose to use a very restricted data set, in order to de- 
rive rather conservative and robust results. 

More generally, we point out that our COSMOMC 
module for computing the primordial spectra nu- 
merically can be used in different contexts, within 
COSMOMC or separately, and even (after minor mod- 
ifications) for studying more complicated models 
producing characteristic features in the primordial 
spectra. The module was written in a user-friendly 
way and made publicly available on the website 
http : //wwwlapp . in2p3 . f r/~lesgourgues/ inflation/ 
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